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q ! ABSTRACT 

We present a detailed analysis of two fields located 26 kpc (~ 5 radial scalelengths) 
<y} | from the centre of M31 along the south-west semimajor axis of the disc. One field 

samples the major axis populations - the Outer Disc field - while the other is offset by 
~18' and samples the warp in the stellar disc - the Warp field. The color- magnitude 
diagrams (CMDs) based on Hubble Space Telescope Advanced Camera for Surveys 
imaging reach old main-sequence turn-offs (~ 12.5 Gyr). We apply the CMD-fitting 
1 technique to the Warp field to reconstruct the star formation history (SFH). We find 

£C) , that after undergoing roughly constant star formation until about 4.5 Gyr ago, there 

CN| ■ was a rapid decline in activity and then a ~1.5 Gyr lull, followed by a strong burst 

\ lasting 1.5 Gyr and responsible for 25% of the total stellar mass in this field. This burst 

appears to be accompanied by a decline in global metallicity which could be a signature 
of the inflow of metal-poor gas. The onset of the burst (~3 Gyr ago) corresponds to 
the last close passage of M31 and M33 as predicted by detailed N-body modelling, 
and may ha ve been triggered b y this event. We reprocess the deep M33 outer disc 
field data of Bar ker et al.l (|201l[ ) in order to compare consistently-derived SFHs. This 
reveals a similar duration burst that is exactly coeval with that seen in the M31 Warp 
■ field, lending further support to the interaction hypothesis. We reliably trace star 

' formation as far back as 12-13 Gyr ago in the outer disc of M31 while the onset of 

star formation occured about 2 Gyr later in M33, with median stellar ages of 7.5 Gyr 
and 4.5 Gyr, respectively. The complex SFHs derived, as well as the smoothly-varying 
age- metallicity relations, suggest that the stellar populations observed in the far outer 
discs of both galaxies have largely formed in situ rather than migrated from smaller 
galactocentric radii. The strong differential reddening affecting the CMD of the Outer 
Disc field prevents derivation of the SFH using the same method. Instead, we quantify 
this reddening and find that the fine-scale distribution of dust precisely follows that of 
the H I gas. This indicates that the outer H I disc of M31 contains a substantial amount 
of dust and therefore suggests significant metal enrichment in these parts, consistent 
with inferences from our CMD analysis. 

Key words: galaxies: individual (M31) - galaxies: interactions - galaxies: stellar 
content - ISM: dust, extinction - Local Group - stars: variables: RR Lyrae 



* Based on observations made with the NASA/ES A Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
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Figure 1. Location of the HST fields (small black squares), over- 
plotted o n the surface densit y of RGB stars from the INT/WFC 
survey of llrwin et al.l (|2005r ) in standard coordinate projection. 
An optical image of the inner disc has been overlaid. The inner 
ellipse has a semimajor axis of 2° (27 kpc) and represents an in- 
clined disc with i = 77? 5 and position angle of 38? 1. The outer 
ellipse, of semimajor axis length 4° (55 kpc), roughly indicates 
the spatial extent of the INT survey. 

1 INTRODUCTION 

Disc-dominated galaxies account for a sizeable fraction of 
the stellar mass in the Universe (e.g., Galla zzTet al.l 12008: 
Ivan der Wei et all 1201 ll ) and thus understanding their for- 
mation and evolution is of prime importance. In the classical 
picture, a disc galaxy forms from the collapse of a rotating 
gaseous protogalaxy within the gravitational potential well 
of a dark matter halo. As the gas cools and settles into 
a rotationally-suppor ted disc, it can fragmen t and undergo 
star formation (e.g., [Fall &; Efstathioul fl980). Modern-day 
studies place these ideas within a cosmological A cold dark 
matter context and use ever more advanced simulations to 
follow the aquisition of mass through mergers with other sys- 
tems, as well as from the smooth accretion of intergalactic 
gas. 

Despite the unprecedented numerical resolution and the 
impressive range of gas physics incorporated, simulations 
have had difficulties successfully reproducing late-type mas- 
sive spiral galaxies such as the Milky Way and M31 (e.g., 
iGovernato et al.ll2009l ; lAgertz, Tevssier.fe Moore 2011) and 
are only now starting to achieve it (fGuedes et al.l l201lf ) . 
A key issue is the so-called 'angular momentum problem' 
which leads to overly small, centrally- concentrated stellar 
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Table 1. Summary of Observations 



Field Warp Outer Disc 



R.A. (J2000) 


00:38:05.1 


00:36:38.3 


Deck (J2000) 


+39:37:55 


+39:43:26 


Rproj^ (kpc) 


25.4 


26.4 


Rdisk^ (kpc) 


31.5 


26.4 


E(B-V£ 


0.054 


0.055 


Dates 


2003 Jun 9-12 


2003 Aug 13-16 


tF606w£ (s) 


8x1310=10480 


8x1250=10000 


tF814U^ (s) 


4x1250+8x1330=15640 


4x1200+8x1280=15040 



a Projected radial distance. 

^ Deprojected radial dista nces are calculated assuming 

(m- M)n=24.47 (lHollandlll998h and an i nclined disc with P.A. = 

38? 1 (iFerguson et al.l2002h and i = 77? 5 (iMa Peng. Gulll997T ). 

c Values from lSchlegel Finkbeiner, fe David (|l998h . 

^ Individual and total exposure times in the F606W band. 

e Individual and total exposure times in the F814W band. 



discs dominated by large bulges. Proposed solutions gener- 
ally invoke various forms and degrees of stellar feedback that 
suppress star formation at early epochs thereby allowing ex- 
tended discs to form, albeit at relatively recent times (z<l). 
Due to the expected inside-out growth of the disc, such mod- 
els require young mean stellar ages at large radii and inter- 
esting constraints arise from observati ons of old and inter- 
mediate age stars in these parts (e.g.. lFe"rguson &; Johnso"nl 
l200ll : IWvsell2008l ). 

Once in place, various processes may influence the 
star formation and chemical enrichment history of galac- 
tic discs. Near or long-range interactions with neighbour- 
ing galaxies can redistribute gas and stars and trigger 
bursts of star formation in which signifi cant fractions of 
the t otal stellar mass can form (e.g., iDi Matteo et al.l 
120081 : IWong et al.|[201ll ). Various authors have also shown 
how stars can undergo large radial excursions due to 
scattering off recurrent transient spirals and other fea- 
tures (e.g., ISellwood fe Binnevl I2OO2I ; iRoskar et ail l2008al : 
iMinchev et al.l l201lh as well as to perturbations cause d 
by satellite accretion (Bird, Kazantzi dis, &; W einberg 201lh . 
The primary effect is to efficiently mix stars over kilopar- 
sec scales within discs, the consequences of which could be 
quite profound. For example, if considerable radial migra- 
tion takes place in discs, it would pose a severe problem for 
interpreting analyses of the resolved fossil record since the 
current location of an old star may be far removed from its 
actual birth place (Ro skar et al.ll2008bh . However, as of yet, 
there is little direct and unambi guous obse rvational evidence 
for radial migration (e.g., [Haywo od 2008). 

The best place to test ideas about disc formation and 
evolution is in the nearest external systems where stel- 
lar populations with a broad range of age, including those 
formed at cosmologically-interesting epochs (z~l-2), can 
be resolved into individual stars. There has been remark- 
able recent progress in this field, thanks largely to the 
ANGST and LCID collaborations which have obtained high- 
quality Hubble Space Telescope (HST) observations of re- 
solved populations in a variety of nearby galaxies (e.g., 
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Figure 2. Colo r-magnitude diagrams f or the Warp (left) and Outer Disc (right) fields, where selected isochrones and a ZAHB from the 
BaSTI library ( Pietrinferni et al. 2004) are overlaid. The error bars show the mean photometric errors as a function of magnitude. In 
each panel, the projected radial distances and the number of stars in the CMDs are indicated. The contour levels correspond to [30, 
55, 80, 105, 130, 155, 180] xlO 3 stars mag -2 . The 50% completeness limit from the artificial stars tests (see Section [274)) is also shown 
by the dashed line in the left panel. The arrows in the top left of each panel indicate the amplitude of reddening (black arrows, from 
ISchleeel. Finkbeiner. fc D avis 1998) and differential reddening (gray arrow, see Section |4j. Note the large spread in the features of the 
Outer Disc field due to strong differential reddening. 



Bernard et aD 120091: iDalcanton et all 120091; IWilliams et all 
2009a; Monelli et al. 2010; Hidalgo et al. 2011). Of particu- 



lar importance are studies of populations in the outermost 
regions of spiral discs: although these regions contain only a 
minor fraction of the total disc mass, they are where many 
evolutionary models diverge the most in their predictions 
and are therefore most easily tested. Interestingly, while the 
star formation histories (SFHs) repo rted thus far show evi- 
dence of inside-out disc growth (e.g., Williams et al. 2009b, 
l20ld : lGogarten et alJl20ld : lBarker et al.ll201lh . thev also in- 
dicate that a significant fraction of the mass in discs was in 
place by z~l-2. 

Over the past few decades, HST color-magnitude di- 
agrams (CMDs) of the outer regions of the our largest 
Local Group neighbours, M31 and M33, have been ob- 
tained and analysed in the context of their SF Hs. Some 
of th ese observations have been limited in depth (Williams 
2002) while others have had their interpretations compli- 
cated by the presence of multiple structura l components 
and/or halo substructure (iBrown et alj|200d: [Barker et al.l 
120071 : iRichardson et al"]l2008l : IRichardson! 1201(1 ). Studies of 
fields sampling the entire M33 disc have revealed a clear de- 
tection of an inv erse age gradient which reverses across th e 
disc truncation ([Williams et al.l [2009b ; Barker et al. I l201lh . 
During Cycle 11, 47 HST orbits were dedicated to observe 
eleven fields in the outskirts of M31 with the Advanced Cam- 
era for Surveys (ACS; P.ID. 9458, P.I.: A. Ferguson). Two 
outer disc fields were observed for ten orbits each, while 
three orbits were secured for each of the remaining halo 



substructure fields. The analysis of the stellar populations 
of the substructure based on this dataset has been presented 
in an earlier series of papers. iFerguson etaD (|2005h explored 
the stellar populations associated with five prominent stel- 
lar overdensities in the halo of M31, and found large-scale 
population inhomogeneities in terms of age and metallic- 
ity. A more in-depth analysis of one of these fields, the 
Gl Clump, placed constraints on the SFH of this region 
and noted the sim ilarity to the M31 outer dis c population 
(|Faria et al.1 l2007h . Finally, IRichardson et al.1 (|2008h com- 
bined this dataset with other deep archival datasets in order 
to explore the nature and origin of the substructure. The 
global picture that emerges from their homogeneous analy- 
sis of 14 HST CMDs is that the substructure can be largely 
explained as either stripped material from the Giant Stellar 
Stream progenitor (|lbata et al.l l200lb or material torn off 
from the thin disc through disruption. 

Here, we focus on the analysis of the two deep disc 
pointings sampling the Warp and Outer Disc. While the 
Warp field was previously presented in IRichardson et al.1 
( 2008), it was analysed at a only fraction of its total depth 
(three out of the ten orbits) and the observations presented 
here are roughly a full magnitude deeper. The paper is struc- 
tured as follows: in Section [2j we describe the observations 
and data reduction, and present the resulting CMDs. The 
CMD-fitting method used to recover the SFH is described in 
Section [3l along with the results for the Warp field in M31 
and the SI field in M33. In Section [H the differential redden- 
ing affecting the CMD of the Outer Disc field is measured, 
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Figure 3. Results from the artificial star tests for the Warp field. 
Top: Completeness in the F606W (green, solid line) and F814W 
(red, dashed line) bands as a function of injected magnitude. 
Middle (bottom) panel: residuals between injected and recovered 
F606W (F814W) magnitudes as a function of injected magnitude. 
The crosses and lines indicate the median and standard deviation 
of the residuals calculated in 1 mag bins. 



and used to constrain the fine-scale distribution of dust in 
the field. The multi-epoch nature of the datasets allowed us 
to search for variable stars in both fields; they are briefly 
described in Section [5] Our interpretation of the results is 
discussed in Section [6] and a summary of the main results is 
given in Section [71 



2 THE DATA 
2.1 Observations 

The two disc fields have been observed with the Ad- 
vanced Camera for Surveys onboard the HST. Their loca- 
tion was chosen so that the on- and off-axis stellar popu- 
lations could be sampled and compared. Throughout this 
paper, we assume a distance of 783 kpc, corresponding 
to (m-M) =24.47 dHollandl Il998h. and an inclined disc 
with PA. 38? 1 feguson et al] l20Q2h and i = 77?5 
(|Ma. Peng, fe Gul Il997h . Located ~18' apart, the 'Outer 
Disc' field samples the outer HI disc along the major axis at 
projected (deprojected) radius of ~ 26 (26) kpc, while the 
'Warp' field samples the stellar disc at a projected (depro- 
jected) radius of ~ 25 (31) kpc, where it strongly bends away 
from the major axis. Note that the deprojected radii, which 
correspond to circular radii within the disc plane, are only 
valid in the case of a planar disc, which is clearly not true 



in the outer regions of M31. Assuming a dust- free exponen- 
tial d isc scalelength of Rd = 5.3 ± 0.5 kpc (| Court eau et al.l 
these radii correspond to ^5-6 scalelengths. 
The fields are shown in Figure [TJ overplotted on a map 
of the s urface density of R GB stars from the INT/WFC 
survey of Irwin et al.l ([2005). For each field, 8 and 12 images 
were secured in the F606W and F814W bands during about 
four consecutive days. Exposure times of ^1,300 seconds per 
exposure lead to a total of about 10,000 and 15,000 seconds 
in F606W and F814W, respectively. A detailed summary of 
observations including coordinates, foreground color excess, 
projected and deprojected radial distances, and exposure 
times is given in Table [1] 



2.2 Photometry and Calibration 

The pipeline processed, non-drizzled ('FLT') data products 
were retrieved from the archive and each exposure was split 
into two individual images corresponding to the two ACS 
chips. Additional processing to correct for intra- field vari- 
ations of the pixel areas on the sky and to flag bad pix- 
els was done using the pixel area maps and data quality 
masks, respectively. Stellar photometry was carried out with 
the standa rd DAOPHOT /ALLSTAR/ALLFRAME suite of 
programs ([Stetson! 1 19941 ) as follows. We performed a first 
source detection at the 3-cr level on the individual images, 
which was used as input for aperture photometry. From this 
catalog, 300 bright, non-saturated stars per image were ini- 
tially selected as potential PSF stars. An automatic rejection 
based on the shape parameters was used to clean the lists, 
followed by a visual inspection of all the stars to remove the 
remaining unreliable stars. We ended up with clean lists con- 
taining at least 200 good PSF stars per image. Modelling of 
the empirical PSF with a radius of 10 pixels (0.5") was done 
iteratively with DAOPHOT: the clean lists were used to re- 
move all the stars from the images except PSF stars, so that 
accurate PSFs could be created from non-crowded stars. At 
each iteration, the PSF was modeled more accurately and 
thus the neighbouring stars removed better. Every 3 to 5 it- 
erations, the degree of PSF variability across the image was 
also increased, from constant to linear, then quadrat ically 
variable. 

The following step consisted of profile-fitting photome- 
try on individual images using ALLSTAR with the empirical 
PSFs previously created. The resulting catalogs were used to 
calculate the geometric transformations between the individ- 
ual images using DAOMATCH/DAOM ASTER. From these 
transformations, it was possible to create a median, mas- 
ter image by stacking all the images from both bands using 
the stand-alone program MONTAGE2. This master image, 
clean of cosmic rays and bad columns and much deeper than 
any individual frame, was used to create the input star list 
for ALLFRAME by performing a second source detection. 

The output of ALLFRAME consists of a catalog of PSF 
photometry for each image. A robust mean magnitude was 
obtained for each star by combining these files with DAO- 
MASTER, keeping only the stars for which the PSF fitting 
converged in at least (N*/2) + l images per band, where N* 
represents the total number of images per band. This con- 
straint greatly limits the number of false detections due to 
noise at faint magnitudes. 
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Figure 4. Color-magnitude diagrams for the Warp and Outer Disc fields showing the location of the regions used for the CMD-fitting 
(left) and for estimating the differential reddening (right, see Section [4|). In the left panel, the dashed boxes indicate the location of the 
two additional bundles used to test our SFH calculation (see Figure [A2]) . 



The final photometry was calibrated t o the VEGA- 
MAG system following the prescriptions of ISirianni et al.l 
(2005), with the revised ACS zeropointfl Note that the 
CMD-fitting technique used in this paper minimizes the im- 
pact of the uncertainties in distance, mean reddening, and 
photometric zero-points on the solutions by shifting the ob- 
served CMD with respect to the artificial CMD; therefore, 
any small errors in these quantities do not affect the results. 

2.3 Color-Magnitude Diagrams 

Figure [2] shows the CMDs for the Warp (left) and Outer Disc 
(right) fields. These were cleaned of non- stellar objects using 
the photometric parameters given by ALLFRAME, namely 
o"606,8i4 ^0.3 and — 0.3 ^ SHAR P ^0.3. Selected isochrones 
from the BaSTI library ([Pietrinferni et al]|2004h are shown, 
both to indicate the location of the oldest main-sequence 
turn-offs (MSTO) and as a guide for compar i ng the two 
CMDs. As already noted by Richardson et al. (2008|), the 
horizontal-branch (HB) of the Warp field is very sparsely 
populated. A theoretical zero-age horizontal-branch (ZAHB) 
for [Fe/H] =— 0.7 from the same stellar evolution library is 
overplotted to sho w its expected location. The redde ning 
(black arrows, from lSchlegel, Finkbeiner, fc Davisl ll998) and 
differential reddening (gray arrow, see Section 3|) at these 
locations are also shown in the top left of each panel. 

Despite their apparent proximity to each other in the 
south-western part of the disc, the CMDs of the two fields 
present notable differences. While the main sequence (MS) 

1 pos ted 2009 May 19 on the HST/A CS webpage: 
:/ /www. stsci.edu/hst/acs/analysis/zeropoints 



of the Warp field contains very few stars younger than ~ 
250 — 300 Myr, the Outer Disc field harbors stars as young 
as a few tens of million years. On the other hand, a plume 
of stars at 0.3 < F606W-F814W ^ 0.6 and F814W > 25.5, 
indicative of strong star formation until about 1 Gyr ago, is 
far more prominent on the CMD of the Warp field. 

Concerning the evolved stages of evolution, both fields 
exhibit well-populated red giant branch (RGB) and red 
clump (RC) features. There are also hints of AGB bumps at 
F814W-23.2. The Outer Disc RGB has a much larger color 
dispersion than that of the Warp, and the fact this spread is 
along the direction of the reddening vector, especially appar- 
ent from the morphology of the red clump (RC), indicates 
that the field is strongly affected by differential reddening. 
We will discuss and quantify this reddening in Section [4] 

Finally, the photometry of the Warp field is about 
0.3 mag shallower than that of the Outer Disc field. This 
is mainly due to the higher background at the time of ob- 
servation. 

2.4 Artificial Star Tests 

In order to retrieve an accurate SFH from a CMD, a good 
knowledge of the external factors affecting the photometry 
is necessary. This is usually done through extensive artifi- 
cial star tests, whereby a large number of artificial stars are 
added to the original images, and the photometry is repeated 
exactly the same way as for the real stars. The comparison 
between the injected and recovered magnitudes reveals the 
biases due to observational effects, such as signal to noise 
limitations, stellar crowding, and blending, as well as quan- 
tifies the completeness of the photometry. 

Unfortunately, these tests cannot quantify the effect of 
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Figure 5. Summary of the 888 calculated solutions in the color- 
magnitude plane: 24 solutions were calculated at each location 
indicated by a plus symbol or a circle. The contours show the 
resulting % 2 distribution in steps of 0.5. Our best-fit solution is 
the average of the 24 solutions indicated by the filled circle. 

differential reddening, as observed in the Outer Disc field 
and described in the previous section. In this particular case, 
the variations are so strong that the CMD-fitting technique 
used to recover the SFH described in Section [3] would be in 
vain. Hence, the artificial star tests were only performed on 
the images of the Warp field. 

The list of stars added to the i mages was produced 
by IAC-star (|A paricio &; Gallartll2004l ). It consists of a syn- 
thetic CMD with very wide ranges of ages (0-15 Gyr) and 
metallicities (Z=0. 0001-0. 02) in order to cover larger ranges 
of color and magnitude than for the observed stars. The 
stars were added on a triangular grid, separated by 20 pix- 
els (twice the PSF radius), to avoid overlap between the 
artificial stars and optimally sample the whole field of view. 
Twenty iterations per image, each with 24,395 artificial 
stars, were necessary to obtain a sufficient sampling of the 
parameter space. In total, almost a million unique artificial 
stars were added to the images for these tests. 

Figure [3] presents the results of the artificial star tests. 
The top panel shows the completeness as a function of 
magnitude for each band. For the F606W and F814W 
bands, the photometry is 90% complete down to 27.29 and 
26.85 mag, respectively, but drops to 50% at F606W=28.75 
and F814W=28.18. As shown in Figured this latter limit 
roughly corresponds to the magnitude of a 12.5 Gyr old 
MSTO. The middle and bottom panels of Figure [3] com- 
pare the injected and recovered magnitudes of the artificial 
stars in each band, as a function of injected magnitudes. 
The offset toward positive residuals at fainter magnitudes is 
a consequence of crowding in these fields, which causes stars 
to be recovered brighter than the input magnitudes. 

3 QUANTITATIVE STAR FORMATION 
HISTORY OF THE WARP FIELD 

3.1 Method 

The SFH of the Warp field in M31 was derived through 
the widely-used CMD-fitting technique, in which the ob- 



served CMD is compared with synthetic CMDs created 
using up-to-date models of stellar evolution. In particu- 
lar, we adopt the I AC method d eveloped for the LCID 
Project (e.g., Monelli et al.l 2010), which uses IAC-star 
([Aparicio Gallartl 120041) to generate s ynthetic popula- 
tions, IAC-pop (lAparicip & ; Hidalgo 2009j)to find the best 
solution, and MinnlAC ([Hidalgo et al.ll201ll ) to produce the 
input files for, and process the output files from, IAC-pop. 

The comparison between the observed and synthetic 
CMDs is done by simply counting the number of stars in 
speci fic areas - referred to as 'bundles' ([Aparicio &; Hidalgo! 
2009) - of the CMDs. These areas are chosen on the basis 
of reliability through considering observational effects (e.g., 
signal-to-noise ratio, stellar density) as well as theoretical 
uncertainties in stellar evolution models. To limit the effects 
of the former, we only consider stars on the CMD brighter 
than the 50% completeness level in both bands. We also ex- 
clude stars on the red giant branch brighter than F814W^25 
(i.e., one magnitude below the RC). The choice to exclude 
the latter is mainly motivated by the confusion between the 
various phases of stellar evolution in this part of the CMD 
- RGB, asymptotic giant branch, RC, red HB, and red su- 
pergiant branch - as well as the uncertainties in the the- 
oret ical models of these advanced evolutionary stages (see , 
e.g. JSkillman et al. 2003; Gallart, Zoccali, &; Aparic iol l2005l : 
I Williams et a l. 2009a). Indeed, most SFH analyses that are 
based on deep MSTO photometry do no t include this part 
of the CMD in the f i tting procedure (e.g ., iBrown et al.| [2006 ; 
iMonelli et ai1l2010l : iHidalgo et alJfeOllh although our tests 
reveal that this does not have a very strong effect on th e 
resultant solution (see Figure lA2l and Bar ker et al"] 1201 lb . 
The bundles used for the Warp CMD are shown as solid 
lines in the left panel of Figure |4] 

The bundles are divided uniformly into small boxes, the 
size of which depends on the density of stars and reliability of 
the stellar evolution models. The size of these boxes ranges 
from 0.06x0.04 mag for the bundles covering the MS to 
0.3x0.5 for the bundle located on the red side of the RGB. 
This yields a total number of ~ 10 3 boxes. Since every box 
carries the same weight, the much larger number of boxes 
on the MS also gives a higher weight to this well-understood 
evolutionary phase in the SFH calculations. 

The number of observed stars, and artificial stars from 
each simple stellar population (SSP), counted in each CMD 
box serves as the only input to IAC-pop. IAC-pop then tries 
to reproduce the observed CMD as a linear combination 
of synthetic CMDs corresponding to SSPs, that is, popula- 
tions with small ranges of age and metallicity. No a priori 
age-met allicity relation, or constraint thereon, is adopted: 
IAC-pop therefore solves for both ages and metallicities si- 
multaneously. The goodness of the coefficients in the lin- 
ear combination is measured through the x 2 merit function, 
which IAC-pop minimizes using a genetic algorithm. These 
coefficients are directly proportional to the star formation 
rate of their corresponding SSPs. 

Age and metallicity are not the only parameters that 
can influence the best-fit solution. Uncertainties in the dis- 
tance, mean reddening, photometric zero-points and model 
systematics, as well as the subjective selection of CMD ar- 
eas and definition of the SSPs, also affect the resulting SFH 
solution. MinnlAC was specifically developed to explore this 
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Figure 6. Best-fit SFH solution for the Warp field, obtained with the BaSTI library. The panels show: (a) the SFR as a function of 
time, normalised to the deprojected area of the ACS field, (b) the age-metallicity relation, (c) the metallicity distribution of the mass 
of stars formed, (d)-(e) the Hess diagrams of the observed and best-fit model CMDs, (f) the residuals, and (g) the significance of the 
residuals in Poisson a. The cumulative mass fraction is shown in red in panel (a). See text for details. 



vast parameter space and combine the best solutions at each 
gridpoint to obtain a robust, stable overall SFH. Its two 
main purposes are: i) to automate the production of the in- 
put files to run a large number of IAC-pop processes in order 
to sample the whole parameter space, and ii) to average all 
the solutions obtained at a given distance and reddening to 
reduce the impact of any specific choices. 

We generated a synthetic CMP conta ining 10 7 stars 
using IAC-star (jAparicio fe Gallartl [2004) with the fol- 
lowing main ingr edients. We select the BaSTI stellar 
evolution library ([Pietrinferni et al . 2004), more specifi- 
cally the scaled-solar models with overshooting, within 
age and metallicity ranges of to 15 Gyr old and 
0.0004^ Z ^0.03 (i.e., -1.7^ [Fe/HK0.18 assuming Z Q 



0.0198: 



K 



roupa 



Grevesse &; No els 1993), respectively. We assume the 
2002h initial mass function (IMF), which is a bro- 



ken power law with exponents of —1.3 for stars with masses 
< 0.5Mq and —2.3 for higher masses. iHidako et all (|201lh 
carried out extensive tests of the effect of different IMF 
slopes for their IC1613 and LGS3 analysis, and found the 
best results were obtained with values compatible with the 
Kroupa IMF. They did note, however, that the resulting 
SFH was rather insensitive to realistic changes of the IMF 
slope. In addition, since our data do not reach the lower 
mass range, a different choice of IMF slope only affects 
the normalizati on and not the rela t ive star formation rates 
(SFRs; see also iBrown et al.ll2006l : iBarker et al.ll200/t ). Fi- 
nally, the binary fraction and mass ratio were set to f=0.5 
and q>0.5, respectively, which is in agreement with G and 
M star surveys of the Mi l ky Way field pop ulation (e.g., 
iDuauennov &; Mavorlll99ll; iReid Gizis 1997). According 
to the tests performed bv lMonelli et al.l (J2010) on five galax- 
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Figure 7. Color-magnitude diagram for the reprocessed M33 field SI from lBarker et all (|201lh , where the deprojected radial distance 
and the number of stars in the CMD are indicated. The contour levels, isochrones, and ZAHB are the same as in the left panel of Figure [2] 
to facilitate comparison. The 50% completeness limit from the artificial stars tests is also shown by the dashed line. Right: same as left 
panel, only showing the location of the four bundles used for the SFH calculations. The dotted lines show the bundles used for M31, 
shifted to the distance and reddening of M33. 



ies, the binary fraction only marginally affects the derived 
SFHs. 

As a first step, the age-met allicity plane was divided in 
19 age bins and 10 met allicity bins, producing the equiva- 
lent of 190 SSPs. The choice of bin sizes (A(age): 0.5 Gyr 
between and 4 Gyr old, 1 Gyr beyond; met allicity bound- 
aries at Z=[0.0004 0.0007 0.0010 0.0020 0.0040 0.0070 0.0100 
0.0150 0.0200 0.0250 0.0300]) was made after testing various 
combinations. However, as stated above, the exact choice of 
these simple populations can affect the solution, so the SFH 
was calculated several times, each time slightly shifting the 
age and met allicity limits of each bin, as well as the location 
of the boxes on the CMDs in which stars are counted. As was 
done for the LCID Project, the age and met allicity bins were 
shifted three times, each time by an amount equal to 30% 
of the bin size, with four different configurations: (1) mov- 
ing the age bin toward increasing age at fixed metallicity, 
(2) moving the metallicity bin toward increasing metallicity 
at fixed age, (3) moving both bins toward increasing values, 
and (4) moving toward decreasing age and increasing metal- 
licity. These twelve different sets of parameterization were 
used twice, shifting the boxes a fraction of their size across 
the CMD. These 24 individual solutions were then averaged 
to obtain a more stable solution. The % 2 of this average so- 
lution is simply the average of the 24 individual reduced % 2 5 
which we denote \ 2 • 

This process is then repeated several times after shift- 
ing the observed CMD with respect to the synthetic CMDs 
in order to account for uncertainties in photometric zero- 
points, distance and mean reddening. Here we shift the CMD 
25 times on a regular 5x5 grid within the following lim- 



its: -0.06 < A(F606W - FSUW) ^ 0.06 and -0.15 < 
AF814W ^ 0.15, as shown in Figure [5] These limits are 
roughly twice the uncertainties in reddening and distance, 
and correspond to AE(B — V) ~ 0.06 and AD ~ 55 kpc. 
Each plus symbol in Figure [5] represents one of the nodes of 
this grid, where 24 individual solutions were calculated and 
averaged. This produced a x 2_ma P where the contours in- 
dicate the approximate location of the absolute minimum 
in the distance-reddening plane. A second series of finer 
shifts around this minimum was applied to refine its loca- 
tion; these are shown as the 12 circles in Figure [5] Therefore, 
888 individual solutions were calculated in total. The filled 
circle indicates the location of the averaged solution with 
the lowest reduced x 2 , i.e., our best-fit solution SFH, with 
X 2 min — 3.08. We stress that we do not consider the po- 
sition of the best x 2 m the distance-reddening map as a 
reliable estimate of distance or reddening since photometric 
zero points, model uncertainties, and other hidden system- 
atics will affect its absolute location. 

The uncertainti es on the SFRs were e stimated following 
the prescriptions of Hidalgo et al ] (|201lh . The total uncer- 
tainties are assumed to be a combination of the Gaussian 
uncertainties gg ~ which include errors on distance, redden- 
ing, and photometric zero points, as well as the effect of 
sampling in the color-magnitude and age-met allicity planes 
- and Poissonian uncertainties ap, due to the effect of statis- 
tical sampling in the observed CMD. Given the large number 
of solutions calculated after varying the input parameters, 
the rms deviation of all the 'good' solutions can be assumed 
to be a reliable proxy of gg- Defining x 2 m in an d as the 

X min 

lowest reduced x 2 an d its standard deviation, we consider 
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Figure 8. Same as Figure [6] but for the reprocessed M33 field SI of Bar ker et al.l (|201ll V 



'good' all the individual solutions with % 2 ^ X 2 m in + cr ^2 ? 

regardless of their location on the x 2_m &P- 

The Poissonian errors are determined by varying the 
number of stars in each box of the observed CMD accord- 
ing to Poisson statistics, with all the other parameters held 
fixed, and recalculating the solution. This process is re- 
peated 20 times, and the rms dispersion of these solutions is 
taken as a p. The total uncertainty is the sum in quadrature 
of gg and ap. 

To assess the significance and reliability of the SFH so- 
lution and to understand the possible systematics, we re- 
peated the calculation of the Warp SFH severa l times either: 
i) us ing the Padova stellar evolution library ([Girardi et al.l 
l2000l ). ii) constraining the range of metallicities for young 
ages (<3.5 Gyr), or hi) including the bright end of the CMD; 
as well as tried to recover the SFH of a series of artificial 
CMDs for which the input SFHs were known. The calcu- 
lations were done strictly following the method described 



in this section. The results proved very reassuring and are 
presented and discussed in the Appendix. 



3.2 Results 

Our best-fit solution SFH is presented in Figure [6] The dif- 
ferent panels show: the evolution of the SFR (a) and metal- 
licity (b) as a function of look-back time; (c) the metal- 
licity distribution of the mass of stars formed; (d)-(e) the 
Hess diagrams of the observed and best-fit model CMDs; 
(f) the residual differences between the two, in the sense 
observed— model; and (g) the absolute residual difference 
normalized by Poisson statistics in each bin. The redshift 
scale shown on panel (a) was const ructed assum i ng th e 
WMAP cosmological parameters from ljarosik et all (|201ll ): 
H 71.0 km s" 1 Mpc" 1 , Q A 0.73 and Q M 0.27. 

The most striking feature of the SFH shown in Panel (a) 
is its complexity; it bears little resemblance to the smooth 
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—0.7 is as expected from the morphology of the RGB 
(see Figure [2}, and in excellent agreement with the stel- 
lar metallicity measured in other outer disc and halo fields 
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Figure 9. Top: H i column den s ity di stribution in the south-west 
disc of M31, from lBraun et al.l (|20Q9h . The location of the Warp 
and Outer Disc fields are shown as the small black and white 
squares, respectively. Bottom: Zoom-in of the Outer Disc field. 
The Hi column density (left), and reddening map derived from 
the spread of the RGB stars along the reddening vector (right) are 
shown (see text). The similarity in morphology of the distributions 
is striking. 



SFHs often invoked in the theoretical models of disc evolu- 
tion. Star formation in the Warp field began early on and 
proceeded at a fairly constant rate for the first ^10 Gyr, 
after which the rate declined rapidly. By that point, roughly 
75% of the stars had been formed. Between 4 and 3 Gyr 
ago, the SFR is consistent with being zero. The follow- 
ing ~1.5 Gyr saw a strong burst of star formation, about 
twice the intensity of the SFR at earlier epochs, during 
which the remaining 25% of the stellar mass formed. For 
comparison, the SFR per unit area of this burst at maxi- 
mum is about half of the current value in the '10 kpc ring' 
(Tab atabaei &; Berkhuii sen 2010), the area with the highest 
current SFR in M31. This burst was followed by a low resid- 
ual SFR until the present time, responsible for the sparse MS 
brighter than F814W ~ 26. The derivation of the SFH using 
the Padova library (see Appendix [Aj) is in excellent overall 
agreement with that reported here. In particular, both so- 
lutions recover a strong burst ^2 — 3 Gyr ago, as well a 
lull in activity ^3 — 4 Gyr ago, and the median age of the 
population, defined as the time at which 50% of the stellar 
mass was in place, is ^7.5 Gyr in both cases. 

The age-metallicity relation (AMR) shown in panel 
(b) is rather well constrained, increasing smoothly from 
[Fe/H] ~ —1.2 at 13 Gyr to slightly above solar metallicity 
2 Gyr ago. From panel (c), the peak metallicity [Fe/H] ~ 



both f rom photometry (e.g., Reitzel, Guhathakurta, & Rich 



Rich et al.1 120041: iFerguson etaD 120051: iBrown et al 
Faria et al.1 120071: iRichardson et al.1 l20Q8h an d spec- 



troscopy (e.g. Jlbata et al.ll2005l : ICollins et al.ll201lh . 

A tantalizing feature in the AMR is the apparent de- 
crease in metallicity from slightly super-solar 2 Gyr ago to 
[Fe/H]~ —0.3 at the present-day. This behaviour is also re- 
covered in the solution calculated with the Padova stellar 
evolution library, although not as strong, but it does not 
appear in any of the tests carried out with mock galaxies 
(see Appendix [Aj . We explored fits in which the metallic- 
ity was constrained to lie within a narrow ra nge at recent 
epoch s but this led to larger residuals (see also lBrown et all 
(2006) who noted a similar behaviour in their M31 fields). 
The decline in metallicity in the most recently formed pop- 
ulations therefore appears genuine although it rests entirely 
on just a few bins. Possible explanations for this behaviour 
include a scenario in which these stars have been accreted 
from a low mas s galaxy (e.g., a Leo A- type dwarf galaxy; 
ICole et alJl2007h or that the stars formed in situ from low- 
metallicity gas which was present in this region ~ 2 Gyr 
ago. 

As shown by the Hess diagrams in panels (d)-(g), the 
model CMD corresponding to this SFH matches very well 
the observed CMD. The only significant residuals in the dif- 
ference and significance plots are due to the RC and red HB 
which, as already discussed, are notoriously difficult features 
to fit with current stellar evolution models. Note, however, 
that the stars of the RC and HB have not been used for the 
calculations of the SFH (see Section |3J"|) . 



3.3 Reprocessing the iBarker et al.l (|201lh M33 
data 

In an attempt to explain the low surface brightness substruc- 
ture r ecently-discovered around M33, iMcConnachie et al.l 
(2009) explored the interaction between M31 and M33 us- 
ing a set of realistic self- consistent N-body models. Their 
favoured orbits bring the systems within ~ 50 kpc of each 
other 2-3 Gyr ago, almost exactly coincident with the ma- 
jor burst of s t ar fo rmation we have found in the Warp field. 
I Barker et al.l (|201lh recently presented the SFH obtained for 
a deep HST/ACS field at the edge of the main disc in M33 
(Rdeproj=9-l kpc, or about 4 scalelengths) which revealed a 
star formation enhancement in the past few Gyr similar to 
that observed in M31. We therefore considered it worthwhile 
to further investigate the nature and timescale of the M33 
burst and examine the evidence that the interaction could 
have triggered simultaneous star formation episodes. In or- 
der to facilitate the most robust comparison, w e opted to re- 
analys e the M33 outer disc data (field SI from lBarker et aT] 
the dataset for which is fully described in that pa- 
per) in the exact same manner as the M31 data presented 
here. This bypasses the difficulties inherent in the different 
choices of algorithms, temporal resolution, IMF, and regions 
of the CMDs to be fit, as well as the representation of the 
best solution itself. We thus carried out the photometry and 
artificial stars tests as described in Section [2] for the M33 
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Figure 10. CMDs corresponding to areas of high (left) and low 
(right) differential internal reddening in the Outer Disc field. Each 
CMD is drawn from ~ 10% of the total area. The same number 
of stars is plotted in both panels. An old isochrone ([Fe/H]=— 0.7 
and 12.5 Gyr old, fro m the BaSTI stellar evolution library; 
IPietrinferni et al1l2004j) is overplotted, assuming E(B -V)=0.055 
from the map of ISchleeel. Finkbeiner. Davisl dl998h . 



field and calculated the SFH following the procedure laid 
out in Section l3J"1 

The CMD is shown in Figure [T) the reanalysis of this 
dataset has led t o a gain of an extra y 0.5 mag in depth over 
that presented in I Barker et al.l (|201ll ). Despite being slightly 
more distant than M31 (A(m - M) ~ 0.2, iBarker et all 
the somewhat longer exposure times and lower stel- 
lar density of the field led to a similar photometric depth 
and 50% completeness level as for the Warp field. In the left 
panel, the same isochrones and ZAHB as in Figure [2] shifted 
to the distance and reddening of M33, hav e been over- 
plotte d for comparison. As already noted by IBarker et al.l 
(|201ll) . the oldest MSTO seen in this field is about one 
magnitude brighter than the 12.5 Gyr old isochrones and 
therefore significantly younger than the age of the Uni- 
verse. We assume a distance of 867 k pc (i.e., (m— M)o=24.69; 
iGalleti, Bellazzini, fe Ferrarol l20Q4f) , inclination of 56°, po- 
sition angle of 23° (|Corbelli &; Schneider! 1997 ), and color 
excess EfB-V)=0.044 as in IBarker et al.l (|201lh . 

In the right panel we show the bundles used to recover 
the SFH as thick solid lines. These are virtually the same as 
the ones used for M31 - shown as the dotted lines, shifted 
to the distance and reddening of M33. The main difference 
is the slight shift of the bundle redward of the RGB to bluer 
colors in order to put tighter constraints on the maximum 
metallicity. 

Our best-fit solution for the outer disc field of M33 is 
shown in Figure [8] For a more meaningful comparison with 
the SFH of M31, panel (a) of both Figure [6] and [8] are nor- 
malised to the deprojected area taking into account the in- 
clination and morphology of each galaxy. The overall agree- 
ment between our solution and that reported in IBarker et al.l 
(2011) is excellent: both solutions indicate a dearth of old 
stars and show that less than half of the total stellar mass 



was in place before redshift 0.5. The chemical enrichment 
laws are also very similar, with a mild metal enrichment 
from [Fe/H] ^ — 1.0 about 10 Gyr ago to roughly —0.5 at 
the present time, and a median [Fe/H] ~ —0.7. We do find 
a slightly higher SFR around ^8-11 Gyr ago than that re- 
ported by IBarker et al.l (|2Q1 lh which may result from the 
increased CMD depth in our reanalysis. This finding is in 
better agreement wit h the presence of a few RR Lyrae stars 
detected in this field ([Barker et al.ll201lh . However, the main 
difference between the solutions results from the finer tem- 
poral grid which we have employed here which allows us to 
resolve a strong and relatively short burst of star formation 
which is exactly coeval to that in the Warp field of M31 and 
has the same duration (^1.5 Gyr). We note that the drop in 
metallicity observed in the SFH of the Warp field in M31 at 
about 2 Gyr ago is also visible in the solution for M33. To- 
gether, these findings reinforce the idea that the enhanced 
SFRs in the outer discs of both galaxies could be related 
and we discuss this possibility in more detail in Section [6] 



4 THE OUTER DISC FIELD 

As shown in Section [2] the Outer Disc field suffers from 
significant differential reddening which, by blurring the fea- 
tures along the reddening vector, affects the morphology of 
its CMD. As we will show below, the dust distribution ap- 
pears highly complex in this region with some stars lying 
both in front and behind the dust clouds. Unfortunately, 
this renders futile attempts to correct for the effects of dif- 
ferential reddening in our data, so the method used in the 
previous section to recover the SFH of the Warp field can not 
be applied here. In particular, differential reddening along 
the RGB mimics the effect of higher metallicity and/or older 
ages by shifting the stars to the red, thus worsening the age- 
metallicity degeneracy already present. 

We focus our analysis of the Outer Disc field on quan- 
tifying the amount of differential internal reddening, under- 
standing its origin and constraining the fine-scale distribu- 
tion of the dust. To estimate the amount of internal redden- 
ing as a function of position in the field, we use RGB stars 
brighter than F814W ~ 26.5, as shown in the right panel 
of Figure |4] since their distribution along an old isochrone 
is roughly perpendicular to the reddening vector. This al- 
lows us to keep contamination from other stellar evolution- 
ary phases to a minimum. We also exclude stars from the top 
two magnitudes of the RGB as their color is very sensitive 
to small changes in metallicity, and age to a lesser extent, 
as well as the part of the RGB which overlaps the RC. 

Given that the blue side of the RGB is well represented 
by a 12.5 Gyr old isochrone with [Fe/H]=— 0.7 (see Figure[2]), 
we assume that all RGB stars should have been observed 
at this location if differential reddening had been negligi- 
ble. Under this assumption, we can estimate the amount 
of extinction along the line of sight to a given RGB star by 
measuring the distance, along the reddening vector, between 
the star and this isochrone, and convert this distance to an 
'effective internal E(B-V)'. Note that this represents the in- 
ternal reddening as it is in addition to the foreground value 
of E(B — V) = 0.055 reported in Table [ J and taken from the 
map of ISchlegel. Finkbeiner. fe David (jl998T ). 
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Figure 11. Distribution of effective internal E(B-V) of the Outer 
Disc field, calculated from the spread of the RGB stars along the 
reddening vector (see text). Note that this reddening is in ad- 
dition to the foreground value of 0.055 quoted earlier. Gaussian 
and double-Gaussian fits to the histograms of the low differential 
reddening area and the whole field-of-view, respectively, are over- 
plotted. The thick black curve is the difference between these two 
Gaussian fits. 

In the bottom right panel of Figure [9] we show the re- 
sulting contour map, where darker tones indicate higher ex- 
tinction. It shows that the distribution of reddening is far 
from homogeneous in this small field-of-view (FOV) and 
delineates clumps and filaments with sizes < a few hun- 
dred parsecs. Interestingly, the distribution of reddening is 
almost an exact m atch to the morp hology of the Hi col- 
umn density map ([Braun et al.l l2QQ9h , kindly supplied by 
R. Braun and shown in the top and bottom left panels 
of Figure [9] By comparing the low resolution H I observa- 
tions of iNewton &; Emerson! (|1977f ) with their photometry, 
ICuillandre et al.l ([200 lh also found that the colors of the 
RGB stars were correlated with the column density of H I 
in their 28 / x28 / outer disc field - which includes our HST 
field. They concluded that the relatively large amount of ex- 
tinction is closely associated with the H I gas, and that the 
outer disc of M31 therefore contains substantial amounts of 
dust. Our result shows that this interpretation is still valid 
on the much finer scales sampled here. 

From our contour map, we can now select the areas 
within the FOV where extinction is lowest and highest. We 
pick stars which fall within two rectangular boxes, each cov- 
ering roughly 10% of the total area of the ACS FOV. Fig- 
ure [10] shows the CMDs of these two fields, each containing 
the same number of stars. The isochrone described above 
is overplotted in each panel, shifted to the distance of M31 
and corrected for foreground reddening. In the left panel, 
the RGB is very wide, almost bimodal, with a blue RGB 
nicely fit by the isochrone and a second, redder RGB. The 
presence of both non- and very reddened stars at the same 
spatial location indicates a complex 3-dimensional spatial 
distribution where non-reddened stars possibly lie in front 
of the dust cloud and the very reddened stars behind. This 



is precisely why we cannot correct the photometry for dif- 
ferential reddening: applying a correction derived from the 
contour map of Figure [9] would simply shift the wide RGB 
to the blue, rather than producing a narrow RGB. 

In the right panel of Figure [TOl the isochrone provides a 
good fit to the stars of the lowest differential reddening area, 
which shows that the differential reddening in this subfield 
is more limited. For simplicity, we assume in the following 
that this CMD is representative of the intrinsic spread in 
color due to a combination of photometric errors and ranges 
of age and metallicity. However, we note that the lower RGB 
here is still not as tight as in the Warp field, and thus some 
differential reddening is still likely to be present if we as- 
sume the intermediate- age and old stellar populations to be 
comparable between the Warp and Outer Disc fields. 

Figure [11] shows the normalised distribution of 'effec- 
tive internal E(B— V)', as calculated above, for the low dif- 
ferential reddening area (gray histogram) and the whole 
Outer Disc field (dash-dotted black histogram). Gaussian 
and double- Gaussian fits to these histograms, respectively, 
are also shown. The difference between the two therefore 
represents the additional contribution of the differential red- 
dening to the spread in the CMD, and is shown as the 
thick black curve in Figure [11] Note that the distribution 
of E(B— V) of the low differential reddening area in Fig- 
ure [TT] is not centred on zero, but is slightly offset by ^0.02. 
While this mi ght indicate that the foreground r eddening 
taken from the Schlegel, Finkbeiner, & David (|l998h map at 
this location is underestimated, it is worth stressing that the 
values of the 'effective internal E(B— V)' are only relative to 
a given isochrone which may not be a perfect representation 
of the true local reddening. Our goal here is to quantify the 
differential internal reddening rather than that due to the 
foreground. The distribution under the thick black curve has 
a median of ^0.16; subtracting the 0.02 offset, we find that 
this additional extinction amounts to ^0.14 (corresponding 
to a total reddening of ~ 0.2) and affects about 40% of all 
the RGB starfl 

While we cannot do a detailed SFH reconstruction for 
the Outer Disc field, we can make an important deduction 
about the recent SFH from inspection of the high and low 
reddening CMDs in Figure [TOl Neither of these CMDs show 
evidence for the prominent 1-2 Gyr turn-off population so 
readily apparent in the Warp CMD hence this region ap- 
parently did not undergo a recent burst. Furthermore, the 
youngest populations (~ few tens of Myr) are largely con- 
fined to the highest reddening regions of the FOV. That such 
stars are present in the Outer Disc field and not at all in the 
Warp is not surprising when considering the gas distribution 
(see Figure [9} which shows that the stellar warp region is 
currently devoid of high-column density gas. 



2 Note that there is no feature in the Schlegel maps at the loca- 
tion of the Outer Disc field. This could be due to the low effective 
resolution of these maps which averages over structures, or the 
fact that outer disc dust is too cold to emit significantly in the 
DIRBE bandpasses. 
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5 VARIABLE STARS 



The candidate variable sta rs were extract ed from the pho- 
tometric catalogs using the lStetsonl (|l996h variability index 
J, which uses the correlation in brightness change in paired 
frames to isolate possible variables from constant stars. This 
process resulted in the selection of ~ 180 and 230 candi- 
dates for the Warp and Outer Disc fields, respectively. Most 
of these turned-out to be false detections, owing their non- 
zero variability index to the combination of sparse sampling 
and large numbers of cosmic rays. Inspection of the indi- 
vidual light curves, and the location of the candidates on 
the stacked images and on the CMDs, allowed us to extract 
the candidate RR Lyrae stars. The other bona fide vari- 
ables, mostly eclipsing binary systems, were discarded due 
to the sparse sampling preventing us from finding their pe- 
riods. Our final catalog of variable stars contains 12 and 13 
RR Lyrae stars from the Warp and Outer Disc fields, re- 
spectively. All of these are new discoveries, since no data of 
sufficient depth or sampling has been obtained before our 
observations of this part of M31. Interestingly, no Cepheids 
were detected, despite the presence of young, massive stars 
in the Outer Disc field. 

We performed the period search through Fourier anal- 
ysis (|Horne &; Baliunasll 19861 ). which was then refined by 
hand to obtain smoother light curves using the sam e inter- 
active program as described in iBernard et all (2009). How- 
ever, the small number of datapoints, as well as an incon- 
venient temporal sampling (all the images in one band ob- 
served consecutively, followed by all the images in the other 
band) significantly complicated the process. For some vari- 
ables, several periods could fold our data and give light 
curves of similar quality. We thus used additional param- 
eters to constrain the most likely period: an amplitude ratio 
Af606W I Af8iaw ~ 1.8, and a reasonable location on both 
the period-amplitude diagram and the CMDs. These con- 
straints were determined from the other RR Lyrae stars for 
which the period search was straightforward. While this may 
not be ideal for obtaining purely objective periods, the pro- 
cess allowed us to obtain reasonably reliable estimates given 
the characteristics of the data. It is important to note that 
some of our final periods could well be aliases of the true 
period. 

Sample light curves for fundamental and first-overtone 
mode pulsators from each fi eld are shown in Fig ure [T2l where 
templates from the set of Lavd en et al . (1999) have been 
overplotted. Given the difficulty of measuring accurate pe- 
riods and mean magnitudes, as well as the incompleteness 
due to the limitations outlined above, we do not attempt 
to analyse the detailed individual and global properties of 
these stars. Their locations and approximate properties are 
given here for future reference. 

The properties of the RR Lyrae stars are summarised 
in Tableland 3 for the Warp and Outer Disc fields, respec- 
tively. The first and second columns give the identification 
number and the type of pulsation, with ab and c designating 
the fundamental (RRafr) and first-overtone (RRc) mode pul- 
sators, respectively. The discrimination between types was 
based on the period and light curve morphology. The next 
two columns contain the equatorial coordinates in J2000.0, 
while columns (5) lists the period in days. Except in the few 





V13 



p=0.56 




V14 



p=0.27 



0.0 



0.5 1.0 
Phase 



Figure 12. Sample light-curves of two RRa6 (top) and two RRc 
(bottom), taken from the Warp (left) and Outer Disc (right) fields, 
illustrating the sparse sampling of our data. Black and grey cor- 
respond to the F606W and F814W magnitudes, respectively. The 
star ID is indicated in the lower left and the period in days in the 
lower right. 



cases where a third digit resulted in a significantly smoother 
light curve (mostly RRc), we only give the periods with two 
significant figures. In the remaining columns, the intensity- 
averaged mean magnitude, as well as the amplitude of the 
pulsations, are listed. 

The detection of RR Lyrae variables this far out in the 
disc shows how the variability studies can complement the 
deep CMD analyses: none of the CMDs presented in this 
paper harbor the obvious features that are traditionally as- 
sociated with truly old stellar populations, such as a well- 
defined faint MSTO or an extended HB. Yet, the character- 
istic properties of these variables demonstrate that at least 
some of the stars at F606W-F814W-0.5 and F814W-25 
belong to an old, extended HB. 



6 DISCUSSION 

The CMD-fitting technique applied to the Warp field reveals 
two main episodes of star formation: an extended period of 
roughly constant star formation over the period from 4.5 
to 13 Gyr ago (corresponding to z> 0.5), and a relatively 
short-lived burst which peaked ~ 2 Gyr ago. In the ~1 Gyr 
which separates these episodes, very little star formation 
took place. 

The differential reddening present in the Outer Disc 
field prevents a similar analysis of its SFH. Although we can- 
not deduce anything about the SFH at early times in this 
field, the lack of a prominent 1-2 Gyr old MSTO feature 
in the CMD indicates that it did not experience the same 
recent burst as seen in the Warp. It is interesting to specu- 
late why the two fields could differ in this respect. Firstly, 
the disc of M31 is strongly perturbed at large radius and 
highly inclined with respect to the line of sight. Thus, even 
though the two pointings appear close in projection on the 
sky they may actually sample very different radii. If the 
SFH over the last few Gyr varies strongly as a function of 
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Table 2. Properties of the Variable Stars in the Warp field 



ID 


Type 


R.A. 


Decl. 


Period 


{F606W) 




{F8UW} 


Af814W 






(J2000) 


(J2000) 


(days) 










V01 


c 


00 38 02.64 


+39 36 30.9 


0.46 


26.03 


0.35 


25.63 


0.39 


V02 


ab 


00 38 10.16 


+39 36 42.8 


0.49 


25.40 


1.25 


25.03 


0.67 


V03 


ab 


00 38 03.65 


+39 37 23.9 


0.47 


24.86 


0.65 


24.46 


0.63 


V04 


ab 


00 38 13.85 


+39 37 37.0 


0.60 


25.79 


0.71 


25.39 


0.29 


V05 


ab 


00 38 13.53 


+39 38 03.9 


0.63 


25.50 


0.69 


24.94 


0.43 


V06 


ab 


00 38 09.02 


+39 38 13.8 


0.70 


25.03 


0.81 


24.52 


0.37 


V07 


ab 


00 38 03.12 


+39 38 20.3 


0.50 


25.12 


1.44 


24.75 


0.71 


V08 


ab 


00 38 03.05 


+39 39 07.5 


0.58 


25.00 


0.79 


24.56 


0.46 


V09 


ab 


00 38 13.77 


+39 39 08.2 


0.53 


24.93 


1.18 


24.44 


0.75 


V10 


ab 


00 38 11.26 


+39 39 13.3 


0.78 


24.77 


0.41 


24.31 


0.22 


Vll 


c 


00 38 11.11 


+39 39 14.2 


0.305 


25.68 


0.27 


25.12 


0.32 


V12 


c 


00 38 13.76 


+39 39 16.8 


0.378 


24.92 


0.58 


24.58 


0.27 



Table 3. Properties of the Variable Stars in the Outer Disc field 



ID 


Type 


R.A. 


Decl. 


Period 


{F606W) 


Af606W 


(FSUW) 


Af814W 






(J2000) 


(J2000) 


(days) 










V13 


ab 


00 36 35.63 


+39 42 17.0 


0.56 


25.88 


1.01 


25.15 


0.94 


V14 


c 


00 36 34.10 


+39 42 27.9 


0.269 


26.34 


0.46 


25.76 


0.18 


V15 


c 


00 36 34.29 


+39 43 02.9 


0.35 


26.25 


0.32 


25.73 


0.33 


V16 


ab 


00 36 43.61 


+39 43 10.0 


0.52 


25.07 


1.18 


24.66 


0.77 


V17 


c 


00 36 42.00 


+39 43 20.4 


0.343 


25.06 


0.40 


24.71 


0.24 


V18 


ab 


00 36 28.25 


+39 43 29.9 


0.59 


24.72 


1.03 


24.36 


0.36 


V19 


c 


00 36 36.18 


+39 43 40.6 


0.35 


25.85 


0.69 


25.30 


0.44 


V20 


c 


00 36 34.25 


+39 43 40.9 


0.405 


25.91 


0.21 


25.12 


0.23 


V21 


c 


00 36 30.60 


+39 43 58.5 


0.38 


25.69 


0.51 


25.12 


0.43 


V22 


ab 


00 36 41.28 


+39 44 01.9 


0.54 


25.39 


0.55 


25.01 


0.39 


V23 


c 


00 36 44.44 


+39 44 41.4 


0.39 


24.97 


0.60 


24.62 


0.35 


V24 


ab 


00 36 42.72 


+39 45 10.2 


0.452 


25.15 


1.17 


24.88 


0.77 


V25 


ab 


00 36 40.05 


+39 45 48.8 


0.51 


24.92 


1.01 


24.54 


0.68 



radius in the M31 outer disc, this could explain the dif- 
ferent behaviours observed. Furthermore, if the two fields 
sample different disc radii then they will also have differing 
levels of contamination from other structural components 
(e.g. thick disc, stellar halo) which could impact the overall 
SFHs. There is some evidence that this might be the case 
from examining published kinematical data for RGB stars 
in the M31 outer disc. Although not exa ctly coincide n t with 
our ACS fields, fields Fl and F2 from llbata et all (|2005T ) 
are in good proximity to the Outer Disc and Warp point- 
ings respectively. While field Fl shows a single velocity peak 
with a velocity lag and moderate dispersion (~ 30 km/s), 
field F2 shows the same component plus a dominant cold 
component (~ 9km/s) which exactly matches the expected 
circular velocity of the disc at that location. This suggests 
that the thin disc may be more dominant in the Warp field 
than in the Outer Disc field, with the thick disc dominating 
in the latter. 



In what follows, we discuss the two main episodes of 
star formation in the Warp field separately, assuming that 
the behaviour observed here is representative of the whole 
outer thin disc. 



6.1 Star Formation at Intermediate and Early 
Epochs 

The star formation rate at lookback times of > 4.5 Gyr ago 
(z> 0.5) is roughly constant and does not show salient fea- 
tures (see Figure [6}. While this behaviour may be genuine, 
it could also (partly) result from the lack of resolution at 
faint magnitudes (hence old ages) which tends to smooth 
the resulting SFH (see the tests on mock galaxies in th e 
Appendix as well as the discussion in iHidalgo et alll201lh . 
One consequence of this smoothing is that some star forma- 
tion is always recovered in the oldest age bins, even although 
none may have taken place. Nonetheless, the presence of RR 
Lyrae stars in the Warp field provides support for a popu- 
lation of stars with age > 10 Gyr. 

If the stars in the Warp field formed in situ, the sig- 
nificant median age of 7.5 Gyr at ^ 5 radial scalelengths 
indicates that M31's outer disc was in place early on. This 
argues against a scenario in which M31 underwent a ma- 
jor merger at intermediate epochs after whi ch the disc re- 
form ed, as has been recently suggested by lHammer et al.1 
(|2010h . Other evidence also supports the idea that the 
outer discs of massi ve spiral galaxies f ormed at high red- 



outer discs ol ma ssive spiral galaxies tormed at high red- 
shift (e.g., iFerguson Joh nson 2001; Bro wn et all [2006: 



IWilliams et~a l. 2009a; Munoz-Mate os et al.ll201lh . On the 
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other hand, the outer disc of M33 l acks a significant frac- 
tion of old stars at 3-4 scalelengths ( Wil liams et al . 2009b; 
iBarker et al.ll201lh . This could reflect the effect of down- 
sizing within the disc population, whereby stars in more 
massive discs tend to have formed earlier and o ver a shorter 
time- span than in less massive systems (e.g., ICowie et al.l 
1996). Alternati vely, it may simp l y high light the uniqueness 
of M33 - indeed iGogarten etafl (|2010h find that the outer 
disc of the similarily low-mass spiral galaxy NGC 300 is 
dominated by old stars. 

Another interesting feature is the fact that the domi- 
nant population which formed at early epochs is relatively 
metal-rich ([Fe/H]> —1.3, Figured), at least in the solution 
obtained with the BaSTI stellar ev olution library. This is 
qualitatively similar to the results of Willia ms et all (|2009aT > 
for the outer disc of M81, which is to date the only other 
massive disc galaxy for which the CMD-fitting technique has 
been applied (although their data do not reach old main- 
sequence turn-offs). The solutions for their field, located at 
five scalelengths from the centre of M81, show a constant 
[M/Hl ~ ~9 - ^ a ^ a ^ a S es - Similarly, our renalysis of the 
I Barker et al.l ([201 lh M33 outer disc field supports the lack 
of stars with [Fe/H] < —1.5, consistent with their original 
findings. Therefore, it seems that even the oldest stars in the 
outer discs of galaxies must have formed from pre-enriched 
gas. 

In the discussion above, we have assumed that the 
stars which we see today in the outer disc of M31 (and 
M33) actually formed there. It is possible, however, that 
these stars formed at much smaller galactocentric dis- 
tances, where the metallicity was higher at early epochs 
due to the higher rate of gas recycling, and migrated 
to their current location due to resona nt scattering with 



trans ie nt spiral arms and bar s fe.g., |Sel lwood & Binne 



3V 

9; 



2002 : iRoskar et al.1 l2008al lbl: ISchonrich fe Binnevl 1200! 
Minchev et al.l 1201 lh or other mechanism s jOuillen et al.l 
120091 ; iBird, Kazantzidis, &; Weinberg! 1201 lh . For example, 
in the models of Rosk ar et al.l (|2008br h roughly 50% of so- 
lar neigbourhood stars (i.e. those lying at ~ 2 — 3 scale- 
lengths) have come from elsewhere, with the bulk originat- 
ing at smaller radii and moving outwards. At extreme radial 
locations such as we have probed here (4 — 6 scalelengths), 
this fraction of migrated stars is likely to be even higher on 
account of the low rate of in situ star formation. Indeed, 
both M31 and M33 exhibit spiral arm patterns which ex- 
tend to these large radii, indicating that migration could 
still be an efficient mechanism for moving stars around. The 
net effect of radial migration on outer disc stellar popula- 
tions is to render the SFH an d AMR flat a nd fe atureless 
in these parts (see Figure 3 of IRoskar etHI (|2008bh ). This 
behaviour is in stark contrast to what we see in the Warp 
and M33 SI fields, where the SFHs display sharp features, at 
least back to >5 Gyr ago, and the AMRs increase smoothly 
with time from early epochs to 2 Gyr ago. This suggests 
only a small role for radial migration in populating the M31 
and M33 outer discs at the radii we have probed. For the 
same reasons, it is also unlikely that a significant fraction of 
the stars in these parts have been directly accreted onto the 
disc plane from satellite galaxies. 



6.2 The 2 Gyr Old Burst 

Strong bursts of star formation are often caused by an ex- 
ternal trigger such as an a ccretion event or an in teraction 
with a nearby galaxy (e.g., [Kennicutt et al.l ll987). Indeed, 
the copious amount of substructure observed in the outer 
regions of M31 strongly suggests interactions have taken 
place (e.g., iFerguson et al.l l2002l : iRichardson et al.1 120081 : 
iMcConnachie et alJl20Q9h . 

The well-kno wn giant stream in the south of M31 
(jlbata et al.ll200lh is the fossil evidence of one such event. It 
has been modeled rather successfully as either the remnant 
of a ~10 9 M(7) sate l lite a ccreted less than a billion years ago 
(e.g., Far dal et al.ll200<t ). or a major merger in which the 
interaction and fusion with M31 may have star ted ~ 9 and 
5.5 Gyr ago, respectively ([Hammer et al.ll2010h . Neither of 
these scenarios can account for a star formation burst that 
occured about 2 Gyr ago in the Warp field. 

Another possibility is that the close passage of M33 
could have triggered star formation on large scales in 
M31. Using realist i c self -consistent N-body simulations, 
IMcConnachie et all (|2009h were able to best reproduce the 
morphology of the outer substructure in M33 with orbits 
that brought the systems within ^50 kpc of each other 2- 
3 Gyr ago and that satisfy their current positions, distances, 
radial velocities, and M33's proper motion. In particular, 
they found that M33 reached pericenter about 2.6 Gyr ago, 
which is in excellent agreement with the start of the star 
formation bursts observed in Figures [6] and [8] This is highly 
suggestive of an interaction between the two systems trig- 
gering inward flows of metal-poor gas which fuelled mod- 
erate bursts of star formation. Such behaviour is consis- 
tent with the predictions of state-of-the-art simulations of 
galaxy interactions and mergers, even though most work to 
date has focused on the inf l ow of gas to the central regions 
(e.g iDi Matteo et al.ll2008l . IrWssier. Chapon. fe Bournaudl 
l2O10h . 

Further development of these N-body models (Dubinski 
et al. in preparation) indicates that the gravitational inter- 
action of M33 should have induced perturbations in M31 
in the form of heating and disruption of the outer disc. 
This is consistent with the presence of warps in both the 
Hi and stellar discs, and the disc-lik e stellar populations 
found in many regions of substructure (|Ferguson et al. 2005; 
iFaria et all 120071 : IRichardson et all l2008h . Such an interac- 
tion could also help explain the large fraction of globular- like 
clust ers younger than 2 Gyr in M31 (jFam de Griis, &; Zhoul 
[2010), the enhanced SF R 2-4 Gyr ago observed in all the 
deep HST fields in M33 (IBarker et alJl2007l : IWilliams etHI 



l2009bl: IBarker et al 



201 ll)" and the warps in both the Hi 
(e.g.. Iwright et al.1 Il97a jPutman et al.1 l2009h and stellar 
discs (|McConnachie et al.ll2010h of M33. To put this inter- 
pretation on firmer footing, additional data are required to 
assess just how widespread the ^2 Gyr burst of star forma- 
tion is in the M31 outer disc and address how the timing 
and intensity of the burst varies with location. 



7 CONCLUSIONS 

We have presented the analysis of very deep, multi-epoch 
HST data for two fields located in the far outer disc of M31 
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at ~ 2° (26 kpc) from the centre, corresponding to 5-6 radial 
scalelengths. 

We apply the CMD-fitting technique to the Warp field, 
and find that a significant fraction of the stellar mass 
(~25%) was produced as the result of a strong and rela- 
tively short-lived burst of star formation abo ut 2 Gyr ago. 
The re - analysis of the M33 outer disc SI field of B arker et al.l 
(Hani) and resulting SFH show a similar burst of SF that is 
exactly coeval with that observed in M31. The burst in both 
galaxies is also accompanied by a decline in global metal- 
licity, which suggests an inflow of metal-poor gas into the 
outer disc. These results further support the inference from 
N-body modelling that M31 and M33 had a close interaction 
about 2-3 Gyr ago that enhanced the outer disc SF at that 
epoch, and may also be responsible for the warps in the H I 
and stellar discs of both galaxies and the large fraction of 
clusters younger than about 2 Gyr old. 

Further back, the SFH was roughly constant over the 
period spanning 4.5-13 Gyr. The presence of stars older than 
about 10 Gyr, while surprising at such a large distance from 
the nucleus, is confirmed by the detection of RR Lyrae stars 
in this field. In addition, the sharp features in the SFH and 
the rather well-defined and monotonically increasing AMR 
favour the in situ formation of these stars, rather than an 
outward migration from birthplaces at much smaller galac- 
tocentric radii. The fact that the oldest populations are still 
fairly metal-rich ([Fe/H]> —1.3) suggests pre-enrichment of 
the outer disc material. 

The second field, the Outer Disc field, is strongly af- 
fected by differential reddening which prevents us from ap- 
plying the CMD-fitting technique. However, the pattern of 
reddening inferred from the color of the RGB stars is highly 
correlated with the Hi distribution in the outer disc, sug- 
gesting these clouds contain substantial dust indicative of 
significant metal enrich ment. This confirms earlier work by 
ICuillandre et al.l {2001) but goes further in demonstrating 
how highly structured the dust distribution is on small 
scales. The presence of significant reddening in these parts 
is consistent with the roughly solar metallicity in the past 
few Gyr uncovered by our SFH analysis of the nearby Warp 
field. 

Finally, from a time-series analysis, we find a dozen 
RR Lyrae stars in each of the two fields. The detection of 
these variables supports the presence of truly old stellar pop- 
ulations in the far outer discs of galaxies despite the lack of 
obvious HBs in the CMDs, and shows how variability studies 
can complement deep SFH analyses. 
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APPENDIX A: TESTING THE SFH 
CALCULATION METHOD 

Here we describe a series of tests aimed at understanding 
how reliable the CMD-fitting technique used in this paper is, 
given the photometric properties of our dataset (i.e., depth, 
crowding), as well as revealing possible biases and limita- 
tions of the method. 

The first test consisted of recovering the SFH of the 
Warp field in the exact same way as described in Section O 
albeit using the stell ar evolution library of the Padova group 
(|Girardi et all l2000h . The resulting SFH is shown in Fig- 
ure [AT] The agreement with the SFH obtained with BaSTI 
(see Figure [6]) is excellent. The main features such as the 
burst centered on 2 Gyr old, the very low star formation be- 
tween 3-4 Gyr and the smooth increase of metallicity from 
the earliest epochs to 2 Gyr ago are virtually the same as 
with the BaSTI library, and the recent decrease in metallic- 
ity is also observed, albeit not as strongly. The main differ- 
ence lies in the overall metallicity, which appears systemati- 
cally lower by roughly 0.2 dex in the former solution at any 
age. The sense and amount o f this offset is exactly the same 
as reported in lBarker et al.1 (|201lh , and is likely due to the 
difference in color of the RG B between the two libraries (e.g., 
iGallart. Zoccali. fe Apariciol I2005T ). We note that the SFH 
obtained bv lBrown et al.l (l20Q6h with the Victoria-R egma 
isochrones ([VandenBerg, Bergbusch, &; Dowlerl [2006) for 
another outer disc field in M31 shows very few stars with 
metallicity lower than about [Fe/H] = — 1, which is more in 
line with the SFH we obtained with BaSTI and agree with 
spectr oscopic estimates of outer disc metal licities in the liter- 
ature (|lbata et al.ll2005l : ICollins et al.ll201lh . Another differ- 
ence between the solutions is that the star formation which 
occurred earlier than 5 Gyr is not as uniform in the Padova 
solution as in the BaSTI one, but shows a mild decrease 
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about 8-9 Gyr ago. Given the larger errorbars in the solu- 
tion obtained with the Padova library, however, we do not 
consider this feature to be very significant. 

For the second test we repeat the fit of the Warp CMD 
with the BaSTI library, but this time including the evolved 
stars - as shown by the dashed bundles in Figure [H The 
solution presented in Figure IA2I is qualitatively similar but 
the fit is significantly worse (x 2 =5.27). This is evidenced 
by the significant residuals in the main-sequence shown in 
panels (f) and (g). We believe the simultaneous fit of all the 
features of the CMD produces worse results due to the larger 
uncertainties in the stellar evolution models of the evolved 
stars. 

The remaining tests involve 'mock galax ies' for which 
CMD s were produced using IAC-star ([Aparicio &; Gallartl 
l2004h and invoked SFHs representative of the typical range 
seen within the Local Group. Observational effects were sim- 
ulated using the results of the artificial stars tests described 
in Section [2.41 Figure [A3] shows the simplest possible SFH, 
that of a constant rate between 14 Gyr ago and the present 
time without chemical enrichment. In panel (a), the input 
SFH is represented by the green lines, while the recovered 
SFH is shown as in the previous figure. In panel (b) the 
two green lines show the upper and lower limits of the input 
metallicity range, between which the stars were uniformly 
distributed in each age bin. Ideally, discrepancies between 
the input and recovered SFH should be representative of 
the effects of the observational effects and age- metallicity 
degeneracy. 

Finally, for the remaining two mock galaxies, we test the 
capacity to recover the SFH of a predominantly old galaxy 
similar to th e transitional dlrr/ dSph of the Local Group 
(e.g., LGS3: iHidalgo et al.l [2oTlh , and of a mostly young 
and intermediate - age galaxy like the dwarf irregular Leo A 
(jCole et al.1 l200/t ). both with a mild chemical enrichment 
over the past 14 Gyr. The results are shown in Figure IA4I 
and [A5l 

Overall, the SFHs are recovered very well in all cases. 
Not only are the residuals in the CMDs very low, but the 
SFRs and metallicities are very close to the input values. In 
addition, the AMR is measured correctly, which highlights 
the ability of the method to solve for both age and metallic- 
ity simultaneously. There are small discrepancies however, 
which appear to be systematic: 

• while the correct mean metallicity at a given age is usu- 
ally recovered, the dispersion of metallicity is often overesti- 
mated, in the sense that stars are recovered outside the range 
delineated by the green lines in panel (b) of Figures IA3I to 
IA51 

• this effect is strongest for the most recent SFH (<1 Gyr) 
since metallicity only has a small effect on the color of the 
bright main sequence and is therefore largely unconstrained; 

• the ages recovered for the oldest populations (>10 Gyr 
old) are slightly overestimated. This may in fact be reveal- 
ing the loss of temporal resolution beyond 7-8 Gyr old due 
to photometric uncertainties at faint magnitudes, as illus- 
trated by the burst of Figure IA4I being recovered with a 
lower amplitude and more spread out in time. 
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Figure Al. Same as Figure [6] but recovered with the Padova stellar evolution library. 
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Figure A2. Same as Figure [6] but using two additional bundles to include the RGB, RC, and HB in the fitting, as shown in Figure 0] 
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Figure A3. Mock galaxy with constant SFR, created and recovered with the BaSTI stellar evolution library. The green solid lines in 
panels (a) and (c) show the input SFH, while in panel (b) they delineate the upper and lower limits of the input metallicity range, 
between which the stars were uniformly distributed in each age bin. The green dashed line in panel (a) is the cumulative fraction of the 
input SFH. 
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Figure A4. Same as Figure [A3] but with a dSph-like SFH. 
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Figure A5. Same as Figure [A3] but with a Leo A-like SFH. 
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